Condensation energy in the spin-fermion model for cuprates 



o 

Oh 



X3 

o 



> 

oo 
m 
a\ 
o 

On 



o 



X 



Ar. Abanov and Andrey V. Chubukov 
Department of Physics, University of Wisconsin, Madison, WI 53706 
(February 1, 2008) 

We use the Scalapino- White relation between the condensation energy and the difference between 
the dynamical structure factor in the normal and the superconducting states to compute the conden- 
sation energy in the spin-fermion model. We show that at strong coupling, the extra low-frequency 
spectral weight associated with the resonance peak in the dynamical structure factor in a supercon- 
ductor is compensated only at energies ~ J which are much larger than the superconducting gap 
A. We argue that in this situation, the condensation energy is large and well accounts for the data 
for cuprates. 
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The understanding of the mechanism of superconduc- 
tivity is an important step towards the general under- 
standing of the physics of cuprates. It has been know 
from the studies of BCS superconductors, that the infor- 
mation about the pairing boson can be extracted from 
measurements of the upper critical field. Specifically, the 
increase of the kinetic energy in a superconductor is over- 
compensated by the decrease of the potential energy as- 
sociated with the feedback effect from superconductivity 
on the bosonic mode which is responsible for pairing [Q . 
The energy difference is called a condensation energy Ec 
and is directly related to the measurable thermodynamic 
critical field by Ec — Vq / (Stt) where Vq is the volume 
of the unit cell. 

Recently, Scalapino and White ||^ applied this reason- 
ing to high Tc superconductors. They argued that if the 
pairing is mediated by spin fluctuations, then the differ- 
ence in the dynamical structure factor S{q, fl) between 
the normal and the superconducting states, integrated 
over frequency and momentum with the weighting factor 
(cos Qx + cos Qy) should be positive and of the same order 
as Ec. This yields a relation 
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where a < 1 is a numerical factor which accounts for 
the fact that the condensation energy is smaller than the 
decrease in the potential energy, and (= 2 for YBCO 
and Bi2212) is the number of layers in the unit cell. 

Neutron scattering experiments in YBCO and Bi2212 
demonstrated 1^,^ that in a superconducting state, 
Sscil,^) possesses a resonance peak at momenta near 
Q = (tt, tt) and at frequencies below 2A where A <C J 
is the maximum of the d— wave gap. The integrated in- 
tensity of the resonance peak yields the r.h.s. of (|l]) con- 
sistent with the data on He- However, it is not clear 
a'priori to which extent the contribution from the reso- 
nance peak measures the actual condensation energy in 
a system. The relevant issue here is whether the spec- 
tral weight of the resonance peak is compensated by the 



depletion of the spectral weight in Ssq{q,oj) at energies 
comparable to A and hence at typical |g— Q| ^ A/J-^ 1, 
or the compensation comes from energies comparable to 
J, i.e., from \q — Q\ of order 1. 

In the first case, the geometrical cos qx + cos qy factor is 
nearly constant for relevant q and can be omitted. Since 
/ (PqdnS{q,Q) ~ S'(5+l)/3 where 5 is the average value 
of the on-site spin (this is the sum rule for spin structure 
factor), the r.h.s of (0) then just measures the difference 
in S between normal and superconducting states. In gen- 
eral, this difference is finite due to a possibility for double 
occupancy for which case d— wave superconductivity fa- 
vors a spin singlet, S = state. However, in cuprates 
double occupancy is energetically unfavorable (Hubbard 
U is large) , and the condensation energy should be small. 

In the second case, however, typical q are far from Q, 
and the momentum dependence of the geometrical factor 
cannot be neglected. In this situation, one can expect 
that the condensation energy is not reduced by the sum 
rule constraint and, generally, is of the same order as the 
unscreened contribution from the resonance peak. 

There are two qualitatively different explanations of 
the resonance peak. One was presented by us and 
is based on strong coupling calculations within the spin- 
fermion model. These calculations extent earlier weak- 
couphng results by others We argued that the reso- 
nance peak is in the particle-hole channel and is related 
to the fact that in a superconductor, the damping of a 
spin fluctuation due to a decay into a particle-hole pair 
is strong only at frequencies above 2A, while below 2A 
it is strongly reduced because of a lack of phase space for 
a decay. By a Kramers-Kronig relation, this reduction 
produces a real part of the spin polarization bubble. At 
low energies, this real part scales as w^, i.e spin collec- 
tive modes in a superconductor behave as propagating 
magnons. This behavior obviously gives rise to a peak 
in Ssc{q, fl) a.t fl = fires Of where ^ is the magnetic 
correlation length. 

Another explanation was presented by Demler and 
Zhang ||] in the context of 5*0(5) theory of supercon- 
ductivity. They conjectured that the peak seen in neu- 
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tron scattering is an antibound state in the spin-triplet, 
particle-particle channel at total momentum Q (tt res- 
onance). Below Tc, particle-particle and particle-hole 
channels are mixed, and the antibound state appears as 
a pole in the spin susceptibility. 

Demler and Zhang recently argued [Q that the mea- 
surement of the condensation energy is a way to distin- 
guish between the two theories. They conjectured on 
general grounds that if the it resonance is the correct ex- 
planation, then it is likely that the compensation of the 
peak spectral weight comes from high energies. Their ar- 
gumentation is that since the antibound state is by itself 
insensitive to Tc, it is energetically favorable for a sys- 
tem to undergo a superconducting transition so that a 
TT— resonance can emerge and increase the condensation 
energy. They also argued that in the spin-fluctuation the- 
ory of the peak, the compensation of the spectral weight 
is confined to a vicinity of 2A, and hence the condensa- 
tion energy is small. 

In the present paper we show that this is not the case. 
We compute S{q, O) in the spin-fermion model and show 
that the compensation of the spectral weight associated 
with the resonance peak in fact comes from high ener- 
gies J or, equivalently, from momenta q far from Q. 
We also argue that at strong coupling, the low-energy 
antibound state in the particle-particle channel does not 
exist because fermionic incoherence washes out the upper 
boundary of fermionic dispersion. 

The point of departure for our consideration is the 
spin-fermion model for cuprates. It describes low-energy 
fermions interacting with their collective spin degrees of 
freedom . Of interest here is the form of the full dy- 
namical spin susceptibility. It has been derived in earlier 
studies [0,|lO|, and we just quote the results. Both in the 
normal and the superconducting state, the full spin sus- 
ceptibility can be written as x~^(q, fi) = Xg ^(q)— nq(f7), 
where xo is the bare susceptibility which is made of 
fermions with energies comparable to bandwidth, and 
nq(f2) is the universal (i.e. cutoff independent) contri- 
bution from low-energy fermions. 

The form of xo is the input for low-energy calculations. 
As before |7l|ic[] , we assume that xa is peaked at Q or near 
Q, and has a simple Ornstein-Zernike form i.e. Xo(q) = 

xoCV(i + (q-Qm; 

The universal contribution to the dynamical suscep- 
tibility involves low-energy fermions and therefore has 
to be computed fully self-consistently within the spin- 
fermion model. Near q = Q, one can neglect q de- 
pendence in n (it yields only a small correction to al- 
ready excising dispersion in Xo(ct)) and restrict with 
IlQ{fl) = Hq. The full susceptibility then has the form 
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In the normal state, TLq is purely imaginary and for any 
coupling strength is almost linear in fl: Hq = i\n\/ujsf 
where ujsf oc The deviations from the linear behav- 
ior result from the corrections to the particle-hole vertex 
which ai ^ = oo are logarithmical in to. However, the 
prefactors are small, and the deviations from linearity 
become relevant only in the extremely tiny region near 
^ = cxD which we will not study here. 

Consider now the superconducting state. Here the 
form of Ho is more complex because spin damping is cut 
below 2 A [^||. By Kramers-Kronig relation, this cut in 
ImUa creates i?eHo which, as we said before, gives rise 
to a resonance peak in x"(Q7 ^) below 2A. 

In general, the spin polarization operator in a su- 
perconductor is a sum of bubbles made of normal and 
anomalous Green's functions. Both fermions in the bub- 
ble has to be near the Fermi surface to satisfy the con- 
straint on energy conservation. For q ~ Q, this restricts 
the momentum integration to the vicinity of hot spots- 
points at the Fermi surface separated by Q. As in 0, we 
consider the situation near optimal doping (when pseu- 
dogap effects are weak) and approximate the gap near 
hot spots by a frequency independent input parameter 
A. Under this approximation, the fermionic self-energy 
and the spin polarization operator Ho are given by a 
set of two coupled equations M : 
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Here R = g/{vF^~^) is a dimensionless parameter which 
governs the strength of the spin-fermion coupling (we use 
the same notations as in |^ - g is the effective spin- 
fermion coupling, vf is the Fermi velocity at a hot spot). 
There are numerous reasons to believe that at and be- 
low optimal doping R ^ 1. To shorten notations, we 
included a bare uj term in G^^ffc, lo) into the self-energy. 

We discuss the solution of (p|JJ) below but first con- 
sider what we actually need to compute. Our goal is 
to check how the extra spectral weight in local Ssd^) 
is redistributed compared to the normal state. For this 
purpose, it is sufficient to compute the integral in (|l|) 
without the geometrical cosgx + cos qy factor and just 
check at which scales the sum rule is recovered. 

Without cos^a; + cos^y, the momentum integration in 
the r.h.s. in (^ can be performed exactly, and at T ^ 0, 
we obtain using S{q, fl) = 2x"iq, f^)/(l - e"''"/^) 
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We absorbed xo^^ factor into the redefinition of Ho. 
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where 



F{fl) = arctan 



arctan 



1 - i?enn 



(6) 



/mllo 

We see that the rate of convergence of the r.h.s. of ^ 
depends on the forms of both -Rello and /mllo above 
the superconducting gap. 

We now obtain these forms from Eq. (^). Quahta- 
tively, the solution of this set has been obtained ear- 
Her . Here we present quantitative results for XIq . 

At i? ^ 1, the normal state self-energy has a Fermi liq- 
uid form T,{u}) cx Z~^{uj + iuj\uj\/ (Aujsf)) with Z oc 
(LOsf/gY^'^ at energies smaller than, and at larger fre- 
quencies crosses over into a non-Fermi liquid, quantum- 
critical regime T,(io) cc exp(i7r/4) lo y/g/\uj\. A simple 
experimentation shows that the solution of depends 
on the ratio between tUsf and the measured superconduct- 
ing gap which is A = AZ if cusf ^ A, and A ^ A^/5 
if LOsf ^ A ^ A. In the first case, at typical frequen- 
cies ~ A the system behaves in the normal state as a 
Fermi-liquid, while in the second case, which is more rel- 
evant to optimally doped and underdoped cuprates [0 , 
fermions with cj ^ A ^ Usf display in the normal state 
the quantum-critical, y/uj behavior. 
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FIG. 1. Schematic representation of -F(Sl) from Eq. 
for (a) weak coupling case iOsf ^ A and (b) strong coupling 
case ujsf ^ A ^ A. Here A and A are measured supercon- 
ducting gaps at weak and strong coupling, respectively, fires 
is the frequency of the neutron resonance peak in a super- 
conductor, and LJsf is a typical spin relaxation frequency in 
the normal state. Observe that in both cases, the frequency 
integral of F{Q) is confined to frequencies, which are much 
larger than the measured superconducting gap. 

We now consider the two cases separately. At uisf ^ 
A, the quasiparticle residue just rcnormalizes the su- 
perconducting gap (A —* A), and the system behavior 
is the same as at weak coupling. In this situation 
ImJIn ~ for fl < 2A, while 1 — Rella changes sign 
at rires = 2A(1 - 0(e-"=^/^). Above 2A, the an- 
alytical form for IIq can be obtained in the limit of 
n > 2A. We found Renqifl) w ttA'^ / {nujsf) and 
ImllQin) ^ {n/LUsf) + i2A^/{nLUsf)) log(r2/A). Substi- 
tuting these results into m) we find after a simple algebra 



that below 2A, F{fl) is negative except for a tiny range 
between fi^es and 2A, while above 2A, F{^l) is positive 
and scales as F{^) oc (1/^^) logJl/A for < Ugf and as 
F{n) cx (l/Sl)^ for > ujsf- This behavior is schemati- 
cally shown in Fig. (|l|a). Splitting the integral in (||) in 
two parts, I = Ii + 12, where the first is the integral over 
frequencies up to twice the measured gap and the second 
is the integral over larger frequencies, and performing in- 
tegration we obtain with the logarithmical accuracy 
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We see that the contribution from low frequencies is neg- 
ative - the vanishing of /mlln overshadows the contribu- 
tion from the resonance peak which for cogf 3> A has an 
exponentially small residue. However, the contribution 
from frequencies above 2A is positive and paramctrically 
larger than negative Ii, such that I ^ Ii + I2 > 0. 

It is essential that although typical frequencies in I2 
are fl ~ Ugf ^ A, still I2 converges at w > Ws/ and 
therefore does not depend on system behavior far away 
from Q. Indeed, at w ~ ujsf, typical q—Q are of the order 
of inverse correlation length (see (Q)). The geometrical 
cos (7a; + cos qy factor is then nearly a constant (= —2) for 
all relevant q and hence a nonzero value of the conden- 
sation energy just refiects the fact that in the absence of 
no double occupancy constraint, the average on-site spin 
in the d— wave superconducting state is smaller than in 
the normal state. Notice that the magnitude of / is much 
smaller than xoA which would be the contribution from 
the resonance peak if its residue was 0(1) 

Consider now the opposite case of ujsf ^ A <C 
A where, we remind, A ~ A'^/g ~ A'^/uigj is the 
measured gap in this limit. Below 2A we still have 
Im Hji = 0, but the resonance frequency (the one at 
which Re Hji ~ 1) is now much smaller than the gap: 

^res ^ (AWs/)l/2 < A l?!. At n > A but Still 

g, we found from ^) that Im Hq approaches 
the normal state form |f2|/ws/, but Re Hq saturates at 
Re Hji = {7TA/{2ujsf){i + 7r~"'^log4) and preserves this 
value as long as the fermionic propagator has a non- Fermi 
liquid, ^/uj form. At very large frequencies -ReHo indeed 
decreases, but the decrease begins only at f2 ~ ^max, 
where either the fermionic propagator recovers Fermi- 
liquid behavior, i.e., bare uj term exceeds ^/uJ contribu- 
tion from the self-energy, or typical q — Q in S{q, fl) 
become 0(1), i.e., lattice effects become relevant. For 
g <C vphp, the recovery of the Fermi-liquid behavior 
comes first, and flmax ~ §■ For g ^ vphp (which in 
Hubbard- model language implies U ^ t), lattice effects 
become relevant first, and flmax ^ (vphp)'^ /g ''^ J- 

Substituting the results for Ho into (j|) we find that 
F{i}) is now positive at frequencies below 2A, except 
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for very low < fi^es ^ 2A. However, F{U) is also 
positive above 2A and, moreover, due to a saturation 
in Rellii, it behaves as F{Q) cx l/fl up to ^max which 
again causes a logarithmical behavior of the frequency 
integral. This behavior of F{U) is schematically shown 
in Fig. (|l]b). Evaluating the integral in with the 
logarithmical accuracy we obtained 

h^fhp r^^^ ^ ^ |£ A/3 log %^ (8) 

where (3=1 + 7r~^log4. We see that Ii is positive, 
i.e., the appearance of the resonance peak below gives 
rise to an extra integrated spectral weight below 2A and 
hence yields a positive contribution to the condensation 
energy. As we discussed before, this extra spectral weight 
should be compensated by a depletion of the spectral 
weight at somewhat higher frequencies. We see however 
that due to non-Fermi liquid behavior of the fermionic 
propagator above 2A, this compensation conies from fre- 
quencies larger than ^max- Moreover, the integrated 
contribution from energies between 2A and ^max is log- 
arithmically larger than the contribution from the res- 
onance peak. Further, for g ^ vphp (the case when 
the no double occupancy constraint is almost exact), 
^max ~ ^sfi"^ and hence typical momenta in xC^i ^) 
are {q - QfS,"^ ^ flmax/uJsf ^ or q - Q = 0(1). 
In this situation, the geometrical cos qx + cos qy factor 
in (|l|) cannot be approximated by a constant and effec- 
tively reduces the contribution from high energies. In 
other words, even if J (PqdilS{q, ft) does not change be- 
tween normal and superconducting states, there is still a 
finite, positive condensation energy which is even larger 
than the net contribution from the resonance peak. This 
is the central result of the paper. 

The magnitude of the contribution to Ec from the res- 
onance peak in optimally doped YBCO has been esti- 
mated in Q without invoking any theory but rather using 
the experimental results for J x"{q,^) They found 
Ec ~ O.OSaJ where, we recall, a < 1 accounts for the 
reduction of E^ due to the increase in the kinetic en- 
ergy. Using J ~ 1500, one obtains Ec ^ 45aK which, as 
Demler and Zhang argued Q agrees with £'c ~ 3 — 12K 
extracted from penetration depth and specific heat mea- 
surements. Our results show that the actual magnitude 
of Ec is even higher due to an extra positive contribu- 
tion from frequencies above 2A. Furthermore, we found 
that this extra contribution to the condensation energy 
is larger than the one from the peak. From this perspec- 
tive, the above estimate for Ec yields a lower boundary 
for the condensation energy. 

To summarize, in this paper we considered the conden- 
sation energy within the spin-ferniion model for cuprates. 
At strong coupling, this model predicts that in a super- 
conducting state, x"{Qi^) possesses a sharp resonance 



peak below twice the maximum of the measured d— wave 
gap. We demonstrated that the appearance of this peak 
does not cause the depletion of the spectral weight in lo- 
cal x" up to frequencies of order J. We computed the 
condensation energy Ec using Scalapino- White relation 
between Ec and x" and found that the dominant, posi- 
tive contribution to Ec comes from a wide range of fre- 
quencies between 2A and J. Our results disagree with 
the assertion in Q that a large condensation energy can- 
not be obtained in the spin-fermion model and therefore 
would require a resonance in the triplet particle-particle 
channel (tt— resonance). 

We on the contrary didn't find any indication of a sharp 
resonance in the tt channel at strong coupling. This reso- 
nance can only emerge as an antibound state and requires 
a sharp upper boundary of the fermionic spectrum. We, 
however, found that strong fermionic self-energy trans- 
forms the spectral weight from the quasiparticle peak to 
higher frequencies and washes out a sharp upper bound- 
ary of fermionic excitations. We caution however, that 
our analysis is valid for a Fermi surface with hot spots. 
Without hot spots, the fermionic decay is forbidden, and 
at least in some range of couplings, the fermionic spec- 
trum preserves a sharp upper boundary in which case the 
system possesses an antibound state in the tt channel. 
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